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Abstract 

The asymptotic properties at future null infinity of the solutions of the 
relativistic Vlasov-Maxwell system whose global existence for small data 
has been established by the author in a previous work are investigated. 
These solutions describe a collisionless plasma isolated from incoming ra- 
diation. It is shown that a non-negative quantity associated to the plasma 
decreases as a consequence of the dissipation of energy in form of outgo- 
ing radiation. This quantity represents the analogue of the Bondi mass in 
general relativity. 

1 Introduction and main results 

The dynamics of a collisionless plasma in interaction with the mean electro- 
magnetic field generated by the charges is described by the relativistic Vlasov- 
Maxwell system. In this model the unknowns are the electromagnetic field 
(E,B) and a set of N non-negative functions f a which give the distributions 
in phase space of N different species of particles. The system consists of the 
Vlasov equation 

d t .fa+p a -d x f a + q a (E+p a AB)-d p f a = 0, Va=l,...,N, 

coupled to the Maxwell equations with charge density p and current density j 
given by 

„ n r N 

p(t,x) = / (YV/ Q )dp, j(t,x)= / (S^q a f a p a )dp. (1.1) 

a =i J ^ 3 a =i 

In the previous equations, i£R is the time, ieR 3 , p€K 3 are the position and 
the momentum of the particles, q a the charge of a particle of species a, 

Pa= -== , p 2 = \p\ 2 
y'm^+p-' 
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denotes the relativistic velocity and m a is the mass of a particle of species a. 
Units are chosen so that the speed of light is equal to unity. 

The relativistic Vlasov-Maxwell system has several applications in plasma 
physics and in astrophysics, where it is used for instance to model the dynamics 
of the solar wind. Many mathematical problems remain unsolved. For example, 
existence of global classical solutions is known only under certain restrictions on 
the size of the initial data (see [21 03 EH global existence for large data has 
been proved for a modified version of the system in which the particle density 
is forced to have compact support in the momentum (see 0]). 

An important feature of the dynamics, which is due to its relativistic char- 
acter, is the presence of radiation fields. The radiation is defined as the part of 
the electromagnetic field which carries energy to null infinity. It is distinguished 
in outgoing radiation, which propagates energy to the future null infinity ^ , 
and incoming radiation, which propagates energy to the past null infinity J^ _ . 
The latter can be interpreted as a flux of energy flowing in onto the system 
from . For an isolated system the incoming radiation should be ruled out 
by appropriate boundary conditions, which will be now briefly discussed. 

The amount of energy £- m (vi,V2) carried to by the incoming radiation in 
the interval [fi,^] of the advanced time, u = £+|x|, can be formally calculated 
by the limit 



where k — x/\x\ and S is the Poynting vector, S— {Ai:)^ 1 {E AB). It should be 
emphasized that (|1.2J) is only a formal definition, since it is not known in general 
whether the above limit exists for a solution of the relativistic Vlasov-Maxwell 
system. Analogously, the limit 



gives the energy which is propagated to ,J? + by the outgoing radiation in the 
interval [1(1,112] of the retarded time, u = t— \x\. 

A solution of the relativistic Vlasov-Maxwell system is isolated from incom- 
ing radiation if £i u (vi,V2) = 0, for all V\,V2 S R- In |2] it was proved that for small 
data this system admits solutions which satisfy this property. These solutions 
are defined by replacing the Maxwell equations with the retarded part of the 
field. The resulting system has been called retarded relativistic Vlasov-Maxwell 
system and reads 



d t fa+p a -d x f a +q a (E rct +p a AB Ic t)-dpf a = 0, Va = l,...,N, (1.4) 




(1.2) 




(1.3) 




(1.5) 




(1.6) 
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where p and j are given by (|l.lfl . The purpose of this paper is to derive infor- 
mation about the asymptotic behaviour at future null infinity of such isolated 
solutions, i.e. to study the properties of the outgoing radiation generated by 
the plasma. Let us first recall, for sake of reference, the global existence result 
of |2. Define 

V = sup{ \p\:(x,p)e supp/ Q (t),teM,l<a< N} 

and denote by -Br(O) the sphere in R 6 with center in the origin and radius R > 
and by A the set of constants {q a ,m a }. 

Theorem 1 Let f a n (x,p)>0 be given in C 2 (R 3 xl 3 ) such that f a n = for 
(x,p,a)eB R (Oyx{l,...,N}. Define A = £ E|p|=o ||V/^||oo, where ,ieN 6 
is a multi-index. Then there exists a positive constant e = e(R,X) such that 
for A < e the retarded relativistic Vlasov- Maxwell system has a unique global 
solution {f a } &(C 2 ) N which satisfies f a (0,x,p) — f a n (x,p) . Moreover F lct — 
(E Tet ,B Tet ) £ C 1 (1R x K 3 ) and there exists a positive constant C = C(R,X) such 
that V < C and the following estimates hold for all (t,x) Gl x M 3 : 

l-FratC*,^)! < CA(1 + 1*| + |a:|)- 1 (H- [t- I^H)- 1 , (1.7) 

p{t,x)<CA(l + \t\ + \x\)- 3 . (1.8) 

The estimate (f 1 . T|l shows that the solution of theorem^is isolated from incoming 
radiation in the sense specified above. We note that the statement of theorem 
n differs from the main result of [2] in two aspects. Firstly in [5] only the case 
of a single species of particles is considered. However the restriction to this case 
has been made only to simplify the notation and the generalization of the result 
of |2] to the case of a mixture is straightforward. Secondly we claim here that 
the distribution functions are twice continuously differentiable, whereas in [2] 
we only proved that they are C 1 . We shall return to this point at the end of 
the introduction. 

We state now the main results of this paper. Let us define 

e(M) = ^(|£rct| 2 + |£rct| 2 )+V / V<+P^f a d P , 

an a "'K 3 

P(t, X) = i (E Iet A B ret ) + V / Pfa dp, 

the local energy and momentum of a solution of (|1.4I) - (|1.6I) . respectively. We 
also set 

M v {u)= I [e-p-k](u+\x\,x)dx (1.9) 
Jr 3 

and note that M v (u) is non-negative. 
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Theorem 2 Let {f a }G(C 2 ) N be a solution of the retarded relativistic Vlasov- 
Maxwell system with data as stated in theorem^ such that F lct E C 1 (R x R 3 ) 
and 

f a (t,x,p) = 0, VxER 3 :\x\>R + a\t\,a = l,...,N, (1.10) 
for some a£ [0,1). Then the limit 

F Tad (u,k)= lim \x\F Ict (u+\x\,x) (1.11) 

\x\ — >+oo 

exists and is attained uniformly in k = x/\x\ E S 2 and uElC, for all Kcl com- 
pact. Moreover the radiation field F rad = (_E rad , _B rad ) satisfies the following al- 
gebraic properties: 

jT^rad t j^r&d q |^rad| |^>rad 

E Iad .k=B I3d -k = 0, 

^rad AB rad). fc= | £ ;rad|2 > 

Theorem 3 Lei {/a} € (C 2 ) N be a solution of the retarded relativistic Vlasov- 
Maxwell system with data as stated in theorem^and such that F lc t E C 1 (R x R 3 ). 
ylssitme the following estimates hold 

(i) V<C 

(ii) \F tA {t,x)\<C{\^t\^x\)-\\ + \t-\x\\)- 1 

(iii) p(t,x)<C(l + \t\ + \x\)- 3 

for some positive constant C and for all (t,x) E R x R 3 . Then Ai v E C 1 and the 
following equation is satisfied: 

^-M y {u) = -—( \E rad {u,k)\ 2 dk. (1.12) 
du An J S 2 

Note that the conclusions of theorems |21 and apply to the solution of 
theorem ^ However the proofs do not require the solution to be small. In 
fact the proofs of theorems [5] and |3] do not make use explicitly of the Vlasov 
equation either. The only tools which enter into play are the continuity equation, 
d t p + y x ■ j = 0, and the local energy conservation law, d t e + X7 X ■ p — 0, which are 
of course satisfied by any "good" matter model. However since the existence of 
global solutions of the retarded relativistic Vlasov-Maxwell system which satisfy 
the assumptions of theorems|5]and|!|]is known by theorem^ we restrict ourselves 
to consider this case. 

Let us now comment the meaning of the results stated above. In theorem 
12 it is claimed that the retarded field generated by the charge distribution 
is asymptotically null and outwardly directed along the future pointing null 
geodesies. (We recall that an electromagnetic field (E,B) is said to be null if 
E-B — and |S| = |B|, cf. |TI], pag. 322.) This result follows essentially from 
\7\. However, for sake of completeness and to help the reader who is not familiar 
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with the formalism used in [7] , we give in section 2 a complete proof of theorem 
13 adapted to our case. The property (|1.10() . which for the solution of theorem 
2] follows from the estimate V < C, is in turn a special case of the assumption 
made in [7J that the matter has to be contained in a timelike world-tube. 

With regard to theorem [3J it shows that the function Ai v plays in the con- 
text of the retarded relativistic Vlasov-Maxwell system the same role as the 
Bondi mass in general relativity (cf. ^ El)- I n fact, A4 V is non- increasing 
and its variation on the interval [1*1,^2] equals the energy dissipated in form 
of outgoing radiation in such interval of the retarded time, as it follows from 
l|1.3fl and theorem [21 Although the Bondi mass loss formula in general rela- 
tivity is extensively studied, it seems the first time that its generalization to 
plasma physics is considered and that l|1.12f> appears in the literature. It should 
be mentioned, however, that a rigorous mathematical derivation of the Bondi 
formula in general relativity is a much more difficult task (cf. 0E1)- 

Let us now deal with the technical point concerning the regularity of the 
solution. We consider for simplicity the system for a single species of particle 
and denote by f rct the unique global solution for small C 2 data. In [2] we 
stated that / rct 6 C 1 ; here we claim that / re t € C 2 . This gain of regularity can 
be easily understood by appealing to the smoothing effect which was pointed 
out in P| . We recall that the solution of the Vlasov equation can be represented 
as f rct = f in (X(0),P(Q)) where X(s),P(s) are the characteristics of (|1.4|) and 
are given by 

X(s) = x + J P(r)dr, 

P(s)=p + J (E let (r, X) + P A B xet (r, X)) dr. 

It was observed in ||| that the time integral of the field evaluated on the char- 
acteristics is one derivative smoother than the field itself provided that X(s) is 
a timelike curve. The latter condition is satisfied by the solution of in virtue 
of the estimate V <C. Hence the characteristics are C 2 and since f ln is also 
given as a C 2 function, then the solution of the Vlasov equation itself is twice 
continuously differentiable. 



2 Algebraic properties of the radiation field 

In this section we prove theorem Let us denote by (f> any of the components 
of the electromagnetic field, i.e. we set <t> — Ei or E>i and define 



F-- 



— {dtji + d Xi p) for the electric field E, 
(d x Aj)i for the magnetic field B. 



Then FeC 1 (Rxl 3 ) and by means of (|rTTU|) . F(t,x) = for \x\ >R+a\t\, The 
retarded field defined by (|1.5|l . (|1.6() has the form 

<Kt,x)= [ F(t-\x-y\,y)- d " 



At,x) \x-y\ 
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where E a (t,x) = {y GR 3 : \y\ < R + a\t— \x — y\\}, which is a compact set for any 
fixed t S R, i£l 3 . We have the following 



Lemma 1 



lim |x|</>(/u+|x|,x) = / F(u + k-y,y)dy, 

Jn a (u) 



Jn a (u) 

uniformly in (u,k) €K. x S 2 , for all K C R compact, where 

«„(«) = {yeM. 3 : \y\ < (l-a)-\R + a\u\)}. 

Proof: The crucial point to prove this lemma is that for i = u+|x| and for |x| 
large , i.e. where we need to evaluate the function </>, the domain of integration 
E a (t,x) is contained in Cl a {u), a compact set whose measure depends only on 
the fixed u. Let p = |x| _1 and define 

2| 12-1—1/2 _ N 



g(p) = [l-2pk-y + p 2 \y\ 2 } 



\x-y\ 



h(p) = p- l [l- pk-y-(l-2 P k-y + p 2 \y\ 2 )^ 2 ] 
= |x| — k-y — |x — y\, forp^O 

and put h(0) = 0, h'(0) = f(fc-y) 2 - ±|y| 2 , so that heC 1 (R) (here the prime 
denotes the derivative with respect to p). Then we have 

lim x <j){u+ |x|,x) = lim / F(u + k-y + h(p),y)g(p)dy. 

Setting G Ui k(p,y) = F(u + k-y + h(p),y)g(p), we have to prove that 

= 0, 



lim 



/ G u . k (p,y)dy- / G Uik (0,y)dy 

JO. a (u) JCla(u) 



1 a (u) JSl a (u) 

uniformly in (u,k) G /C x S 2 . By the mean value theorem we have 

G u ,h(p,y) = G Utk (0,y)+pRu,k(P>y)> 

where the remainder is bounded as 

\R»,k(j>,v)\ < su P { | G' u ^(T,y) |, < r< p}. 

Hence 

/ G u . k (p,y)dy~- / G«,fc(0,y)dy 
< p sup { \G' uM (t, y) | , < r < p, y e o («) } Vol[ft a (u)] . 
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Since G is C 1 and fl a (u) is compact, the lemma is proved. □ 
By means of lemma 1, the radiation field is continuous and is given by 

E™ d (u,k) = - f (d iP +d t ji)(u + k-y,y)dy, (2.1) 

B™ d (u,k)= [ (d x Aj) t (u + k-y,y)dy. (2.2) 

It remains to prove the algebraic properties of (E r&d ,B rad ). In we replace 
the identities 

dtji {u + k-y,y) = d u j l {u + k-y,y), 

dip(u + k-y,y) = dy t [p{u + k-y,y)]+d y -[3(u + k-y,y)}ki 
-d u (j ■k)(u + k-y,y)k l , 

the second one being a consequence of the continuity equation. After integrating 
by parts we get 

E™ d {u,k)= f d u [(j-k)h~j t }(u + k-y,y)dy. 
Jn a (u) 

Let M denote the vector 

M(u,k)= [ [{j-k)k-j](u + k-y,y)dy. (2.3) 
Jn a (u) 

Since the integrand function in l|2.3|) is C 1 and vanishes on the boundary of 
n a (it), then we have E rad = d u M. Analogously, from $X?fy we get B rad = d u N, 
where 

N{u,k)= [ {jAk)(u+k-y,y)dy. (2.4) 
Jn a (u) 

The next lemma describes the algebraic properties of the vectors M, N. 
Lemma 2 V(«,i;)6Rx5 2 the vector fields defined by \2. \2.J$ satisfy 

(1) M(u,k)-k = N(u,k)-k = 

(2) M(u,k)-N(u,k) = 

(3) \M{u,k)\ = \N{u,k)\ 

U) (M(u,k)AN(u,k))-k=\M{u,k)\ 2 . 

Proof: The proof of (1) is straightforward. For (2) we put j = j(u + k-y,y) and 
j' = j(u + k-y' ,y') for short and write 

M-N = j (j -k)k\- {]' A k) dy' dy 

j-(f Ak)dy'dy 
j'-(jAk)dy'dy, 
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where it is understood that the integrals are over the set Q a (u). Then inter- 
changing y and y' we get M ■ N = —M ■ N, i.e. M ■ N = 0. To prove (3) we write 



\M\ 2 = I I (j-(j-k)k)-(j'-(f-k)k)dy'dy 
\j-j'-(j-k)(j'-k)}dy'dy 



and 

|Af = / / (jAk)(fAk)dy'dy. 

In the previous equation we use the following rule of vector calculus 

(aAb)-(cAd) = (a-c)(b-d)-(a-d)(b-c), 

which is valid for any vectors a, b, c, d and the identity (3) follows at once. To 
prove (4) we write 



(MAN)-k = I \[{{j-k)k-j)A{j' Ak)]-kdy'dy 
[{{j-k)k-j)Ak]-(f Ak)dy'dy 
(jAk)(fAk)dy'dy=\N\ 2 =\M\ 2 . 



□ 

The following lemma permits to relate the algebraic properties of the vectors 
M, N to the ones of the radiation held and concludes the proof of theorem 
To simplify the notation we suppress the dependence on k and denote by an 
upper dot the differentiation with respect to u. 

Lemma 3 Let M(u),N(u) be C 1 vector fields satisfying the properties (l)-(4) 
of lemma 2. Then Vw G K: 

(a) N(u)-k = M(u)-k = 

(b) M(u)-N(u) = 
(e) |M(«)| = |JV(«)| 

(d) (M(u)AN(u))-k=\M(u)\ 2 . 

Proof: The identity (a) is proved at once by differentiating (1) with respect to 
u. To prove the other identities we consider a coordinate system in which the 
z— axis is parallel to fc. In this frame the vectors M and N have the form 

M(u) = (mi(u),m 2 (u),0), 
N(u) = {n 1 (u),n 2 (u),Q). 
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Now, because of (2), (3) and (4) of lemma 2: 

m\{u)ni{u) + 77i2 (it)n2 (it) = 0, (2-5) 

mi (it) 2 +7772 (u) 2 — ni(u) 2 + n 2 (u) 2 , (2.6) 

777l(Tt)772(tt) — 777-2 (ti)t7i (it) = 777l(u) 2 + 7772 ( 7i) 2 . (2-7) 

After some elementary algebra, <|2.5[) - (|2.7|) give 77ii =712, 7772 = —771. Hence the 
vectors M,N can be represented in the following form: 

M(u) = (q(u),r(u),0) => M = (q,r,0), 
N(u) = (-r(u),g(u),0) (-r,g,0), 

by which the properties (b), (c) and (d) follow at once. □ 



3 Bondi mass of the plasma 

In this section we prove theorem Let us introduce 

m v (r,u)= / [e — p-k](u+\x\,x)dx. (3.1) 

J I x I < r 

Note that m v (-,Tt) is non-decreasing and so its limit as exists. We first 

prove l|1.12f> assuming that 

(jj) m v (r,w) converges as for all it£R. 

Assume (jj) holds. Then A4 v (u) is well defined as improper integral and we have 

M v (u)= lim m v (r,7t). 

r — >+oo 

Let JC be a generic compact subset of R. Evaluating the energy conservation 
law dte = —d x -p on the future light cone corresponding to the value 71 of the 
retarded time we have 

d M e(u+ \ x\,x) = —d x -p(u+\x\,x). 

In the previous equation we use the identity 

d x ■ p (u + \x I , x) = d x ■ [{p {u + I x I , x)} - d u (p ■ k) (u + \x \ , x) 

and so doing we get 

d u (e-p-k)(u+\x\,x) = -d x -[p(u+\x\,x)}. (3.2) 

We now integrate l|3.2|l on the region |x| <r, and use the Gauss theorem to 
transform the right hand side into a surface integral over the sphere of radius r. 



9 



Since f(t,x,p) is supported on the region \x\ <R + a\t\, with a£ [0,1) then for r 
large enough, f(u+ \x\,x,p) vanishes on S r — {x : |a;| = r} and so we get 

lim / d u (e — p ■ k)(u+ \x\, x)dx 

r ^+ oc J\x\<r 

= -(Att)- 1 lim / (E Iet AB iet )-k(u + r,x)dS r . (3.3) 



/' — ' + OC 



The integral in the left hand side of l|3.3|l is equal to 9 u m v (r,u), W>0. By 
theorem |21 and l|3.3|l , d u m converges to the right hand side of (I1.12|) . uniformly 
in u G /C, as r — > +oo. Hence, using (jj) we infer that m v (r,w) converges uniformly 
in u € K, and also 

lim d u m v (r,u) — -f- ( lim m v (r,u)) = -^-M^ (u). 

r^ + oo dll Vr^+oo / du 

Thus the proof of (|1.12|) is complete if we show that the property (jj) above 
is satisfied. We note that to this purpose, a direct use of the estimate (ii) in 
theorem |3] is not enough, since it entails \E(u+\x\,x)\—0(\x\~ 1 ), as |a:| — »oo, 
which is too weak to bound A4 y . To solve this problem we rewrite A4 V in a 
way that the decay at past null infinity, which is faster by means of the absence 
of incoming radiation, enters into the estimates (in other words, we make the 
advanced time v appear instead of the retarded time u). For this purpose we 
introduce, besides m v (r,u), the following function 

m A (r,v)= (e+p-k)(v- \x\,x)dx. (3.4) 

We also set 

M A M= lim m A (r,v). (3.5) 

r — >+oo 

Lemma 4 Under the assumptions of theorem^ A4 A (v) is continuously differ- 
entiable and satisfies 

j!-M A (v) = 0, VveM. (3.6) 
dv 

Proof: We split Ai A (v) into four parts as follows: 

M A (v) = h (v) + h (v) + h {v) + h(v), 



where 



h(v) = ^~ [ (l^ctl' + l^etl 2 )^-^!^)^, 

h(v) = y2 / / V m a+P 2 f<x(v- \x\,x,p)dpdx, 

- JR3J\p\<C 

I 3(v) = ~r / (E Ict AB 1 . ct )-k(v-\x\,x)dx, 

47T ./ B3 
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h(v)=y~) I I p-kf a (v~\x\,x,p)dpdx. 



p\<c 



Since I3 (resp. I4) is dominated by I\ (resp. I2), it suffices to estimate h(v) and 
l2(v). By means of (ii), I\{v) is uniformly bounded in ws/C. For ^(w) we use 
that f a (v— \x\,x,p) = for |x| >iZ+o|v — |aj||, which implies / Q (i> — \x\,x,p) = 
for \x\>(l-a)- 1 (R + a\v\). Thus 

|/ 4 (w)|<CVol[{a;:|x|<(l-a)- 1 (i? + a|?;|)}]<C, Vwe/C. 

Since m^(-,v) is non-decreasing, the limit (|3.5|) exists for all d£| and by the 
previous estimates, A4/\(v) converges, as improper integral, uniformly in veJC. 
To prove (|3.6(l we repeat the argument which led to (|3.2() . now evaluating on 
the past light cone. So doing we get 

d v (e+p ■ k){v-\x\,x) = -d x ■ [p(v- \x\,x)]. (3.7) 

Integrating and using the Gauss theorem we obtain, 

lim / d v (e + p-k)(v—\x\,x)dx 



r — s--(-oo 



|a:|<r 



= -(4tt)- 1 lim / (E Iet AB Tet )-k(v-r,x)dS r , (3.8) 

where again we used that f(v — r,x,p) vanishes on S r for large r. The estimate 
(ii) implies that the right hand side of Ij3.8[l tends to zero uniformly in v S K 
and so we get 

0= lim / d v (e + p-k)(v— \x\,x)dx 



r — >+oo 



|x|<r 



d f d 

= — lim / (e+p-k)(v-\x\,x) = —M A (v), 

dVr->+ooJ\ x \< r dv 

where the uniform convergence has been used to shift the derivative. □ 

Note that in the proof of lemma 0] the estimate (iii) in theorem [3] has not been 
used. This will become important for the completion of the proof of (Jj). The 
identity 1)3. (jf) represents the counterpart of l|1.12f) on the backward cones of 
light. In fact, from one hand the conservation of M.^(v) is due to the absence 
of incoming radiation and, on the other hand, the Bondi mass M. y {u) decreases 
as a consequence of the emission of outgoing radiation. 

We are able now to complete the proof of (jj). Integrating (|3.7|) between u 
and u + 2|a;| we get 

ru+2\x\ 

(e+p-k)(u+\x\,x)-(e+p-k)(u-\x\,x) = - / d x -[p(v-\x\,x)]dv. 
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Integrating in the region \x\ < r we get 

(e + p-k)(u+\x\,x)dx — / (e + p-k)(u-\x\,x)dx (3.9) 

x\ <r J |a;| <r 

/•u+2|a| 

/ d x -[p(v—\x\,x)]dvdx. 

\x\<r J u 

Now we use the identity 

®x [ I p(v—\x\,x)dv] — 2p-k(u+\x\,x) 

J u 

pu+2\x\ 

+ / d x -[p(v-\x\,x)]dv, 

J u 

Substituting into l|3.9fl and using the Gauss theorem we obtain 

m v (r,u) —m A (r,u) + q(r,u). (3.10) 

where 

q(r,u) = — I / p-k(v—r,x)dvdS r 

JS r Ju 

1 f r u + 2r 
= — — / / (E let AB let )-k(v-r,x)dvdS r 

u+2r 



r ru-t-zr p 

y~) I I / p-kf a (v-r,x)dvdS r 

„ JS r Ju J\v\<C 



ISrJu J\p\<C 

Using (ii) in the first term and (iii) in the second term, we conclude that 
lim. r ^ +00 |q(r,w)| is bounded for all mg! and so the property (jj) follows from 
H3.10[l and lemma 01 This concludes the proof of theorem 
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